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Abstract 

The lower moments of the unpolarized and polarized deep-inelastic structure 
functions of the nucleon are calculated on the lattice. The calculation is done 
with Wilson fermions and for three values of the hopping parameter k, so that we 
can perform the extrapolation to the chiral limit. Particular emphasis is put on 
the renormalization of lattice operators. The renormalization constants, which 
lead us from lattice to continuum operators, are computed perturbatively to one 
loop order as well as non-perturbatively. 
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1 Introduction 



For our theoretical understanding of the short-distance structure of the nucleon, as 
well as for a successful explanation of the recent HERA and polarized lepton-nucleon 
scattering data, a calculation of the deep-inelastic nucleon structure functions from 
first principles is needed. The theoretical basis for such a calculation is the operator 
product expansion. The quantities of primary interest are the lower moments of the 
quark and gluon distribution functions and their higher-twist counterparts. 

We have initiated a program to compute the moments of the unpolarized, Fi and 
and polarized nucleon structure functions, gi and (/2, on the lattice. First results of 
our calculation have been reported in Ref. [1]. In this talk we shall focus on two topics: 
the valence quark distribution and the spin content of the nucleon. The calculation 
of the gluon distribution functions and the distribution functions involving sea quarks 
is in progress. Because of space limitations we shall not be able to discuss the higher 
moments of gi and the structure function g2. 

Thus we shall concentrate on the moments (x""^), where x is the fraction of the 
nucleon momentum that is carried by the quarks, and Aq^ the quark spin contribution 
to the nucleon spin. We have 

where 

\Y.iP^^\^^{l-^n}\P^^ = 24^^bM ■■■Pf^n- traces], 

S 

^ii^-^n = ( 2 j D ^^■■■D ^„q- traces (2) 

with q = u[d) for / = u[d). Here {• • •} means symmetrization of the indices. Further- 
more we have 

A? = ^ai^\ (3) 

where 

O'J^^ = q^.l.q. (4) 

The lattice calculation is performed on a 16"^ X 32 lattice at gauge coupling fi = 
6/(/^ = 6.0. We work in the quenched approximation where one neglects the effect 
of virtual quark loops. We use Wilson fermions, and we compute simultaneously at 
three values of the hopping parameter, n = 0.155,0.153 and 0.1515, so that we can 
extrapolate our results to the chiral limit. This translates into physical quark masses 
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niq of roughly 70, 130 and 190 MeV, respectively. So far we have collected of the order 
1000, 600 and 400 independent configurations at the three k, values. 

The lattice operators O are obtained from the operators in the euclidean continuum, 
up to factors of i, by replacing the covariant derivative in (2) by the lattice covariant 
derivative ^ 

D ^{x,y) = -[Ui,{x)Sy^^+f, - Ujix - fi)Sy,^-f,]. (5) 
We first compute the two- and three-point correlation functions 

Cr{tJ) = Y.^pABa{t.p)Bp{Q,p)), 

a,fi 
a,f3 

As our basic proton operator we use (with C = 7472 in our representation) 

Ba{t,p)= E e-'P^'e,,,ul{x){u'{x)Cjsd^xj). (7) 

x,a,b,c 

The nucleoli matrix elements of interest are then obtained from the ratio 

= i El — -TrlTNJN] (8) 

(for ^ T ^ t), where N = {E^f4 — ip^ -\- m]\j)l Ej^ and J is defined by 

(p, s\0\p, s) = u{p, s)J u{p, s). (9) 

In order to increase the overlap with the ground state and to make the plateau 
region - by this we mean the region where the excited states have died out - in r as 
broad as possible, we use 'Jacobi smearing' [2]. We smear both source and sink. In 
Fig. 1 we show the effective nucleon mass for n = 0.155, i.e. our lightest quark mass, 
as given by ln(C(t)/C(t -|- 1)). We find a good plateau. Our results for the hadron 
masses are compiled in Table 1. The chiral limit is obtained by extrapolating in Ij n 
to = 0. If we assume that vanishes linearly with 1/k — )• 1/kc, we obtain the 
critical value Hc = 0.15693(4). 
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0.1515 


n 

0.153 


0.155 


nip 
rriN 


0.504(2) 
0.570(2) 
0.900(5) 


0.422(2) 
0.507(2) 
0.798(5) 


0.297(2) 
0.422(2) 
0.658(5) 



Table 1: The hadron masses in lattice units at /3 = 6.0. 
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Figure 1: Effective nucleon mass plot. Both source and sink are smeared. The hori- 
zontal line indicates the result of the fit as well as the fit interval. 
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Observable 


(C) 


Components 


10 


no 


{x) 
Aq 


V2 
V3 
V4 

ao 


0{44} - UO{ll} + 0{22} + 0{33}) 
C{114} — |(C{224} + C'{334}) 
C{1144} + C'{2233} — C'{ii33} — C{2244} 
01 


16 

3 
25 

3 
157 
15 




-1.892(6) 
-19.572(10) 
-37.16(30) 

15.795(3) 



Table 2: The lattice operators. The momentum is taken to be p = (27r/16,0,0) = 
(pi, 0, 0) in the case of V3 and and p = elsewhere. 

In our calculation of the three-point functions we have fixed t at 13. For the 
unpolarized case we have taken V = |(1 -|- 74). For the polarized case we have chosen 

r = 1(1 -1-74)17572, which corresponds to polarization -| in the 2-direction. For the 

calculation of the higher moments we need non- vanishing nucleon momenta. We have 
taken p = and p = (27r/16,0,0) = (pi,0,0). 

2 Renormalization of Lattice Operators 

The bare lattice operators, (9(a), which are in general divergent, must be renormalized 
appropriately before we can use them. We define finite operators (9(/i), renormalized 
at the scale by 

0{^i) = Zo{{a^i)\g{a))0{a), (10) 

where we define 

{q{p)\0{lJ:)\q{p)) = {q{p)\0{a)\q{p)) I^^^V, (11) 

\q{p)) being a quark state of momentum p. In the limit a — )• this definition amounts 
to the continuum, momentum subtraction renormalization scheme. The numbers that 
we will quote later on will all refer to this scheme. 

The lattice operators must be constructed such that they belong to a definite ir- 
reducible representation of the hypercubic group -ff(4) [3, 4]. In particular they must 
not mix with lower-dimensional operators. In this talk we will consider the operators 
listed in Table 2. 

We have computed the renormalization constants of these operators (and others) 
in perturbation theory to one loop order [1, 5]. (See also Ref. [6].) We write 

Zo{{a^i)\g) = 1 - -|^C^[7o ln(«/«) + ^o], (12) 
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Figure 2: The non-perturbatively determined renormalization constant Z^^ as a func- 
tion of . The sohd hne represents the perturbative result given in Table 2. 
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Figure 3: The non-perturbatively determined renormalization constant Z^^ as a func- 
tion of /i^. The solid line represents the perturbative result given in Table 2. The 
dashed line is the result of tadpole improved perturbation theory. 
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where Cf = 4/3 and 70 is the anomalous dimension that also enters the Wilson 
coefficients. The results of our calculation are listed in Table 2. In case of V3 there 
is a small mixing problem [5], which we have ignored here, because numerically it is 
insignificant. 

The effect of renormalization can be relatively large, in particular for those operators 
that involve higher powers of covariant derivatives. The main sources of contribution 
are the leg and operator tadpole diagrams. There has been a lot of discussion in the 
literature [7] on how to reorganize perturbation theory around these contributions in 
order to achieve a better convergence of the perturbative series. 

It is important to determine the renormalization constants accurately. In view of 
this we have computed the renormalization constants non-perturbatively as well [8], 
following the suggestion of Ref. [9]. The renormalization constants are obtained from 
the calculation of the operator matrix element 

{qip)\Oia)\qip))\^2^^2. (13) 

For this calculation we need to fix the gauge. We have chosen the Landau gauge, and 
we average over all Gribov copies. In Figs. 2 and 3 we show our results for Zy^ and Zag 
as a function of jj,^ (in lattice units). Here k = 0.153, but we do not see any dependence 
on K. We compare the non-perturbative results with the perturbative calculation. (The 
renormalization scheme used in the non-perturbative calculation differs slightly from 
the perturbative one. But the difference is insignificant [8].) We find remarkably good 
agreement between the two approaches. In Fig. 3 we have also shown the prediction of 
tadpole improved perturbation theory [7]. The non-perturbative result lies in between 
the perturbative and the tadpole improved result, so nothing seems to be gained by 
tadpole resummation in this case. But we will have to wait until we have computed all 
renormalization constants non-perturbatively before we can draw any firm conclusions. 
In the following we shall take 

= ^ 2GeV', (14) 

which eliminates the logarithms in the renormalization constants, and we will denote 
Zo{l,g = 1.0) hy Zo. 



3 A Selection of Results 

We are now ready to compute the nucleon matrix elements (2) and (4). The physical 
matrix elements at the scale fi are obtained from the ratios (8) by 

(15) 



Zy^ 2k Zao 2k 2Lp^ 
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Figure 4: The moments {x^~^) as a function of together with a hnear fit to the 
data. The sohd circles indicate the extrapolation to the chiral limit. The diamonds 
mark the phenomenological valence quark distribution of Ref. [10] (fit -D-). 
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Figure 5: The quark spin contribution to the proton spin as a function of together 
with a hnear fit to the data. The sohd circles indicate the extrapolation to the chiral 
limit. 

We have made sure that we are computing the matrix elements of the lowest-lying 
state, i.e. the nucleon. We find that the signal is practically constant for distances of 
T larger than two lattice spacings from the source [t = 0) and from the sink [t = 13). 
We have defined the continuum quark fields by \/2k times the lattice quark fields. For 
the renormalization constants we take the perturbative values given in Table 2. The 
results are plotted in Figs. 4 and 5. All our results are given for the proton. The 
distribution functions of the neutron are obtained by interchanging u and d. 

The first observation we make is that the nucleon matrix elements show roughly a 
linear behavior in 

+ l-Sm^a, (16) 

i.e. in the quark mass. The lines drawn in Figs. 4 and 5 are linear fits to the data. 
The result of the extrapolation is indicated by the solid circles. 

Let us next discuss the moments {x^~^). We find that the lowest moment is practi- 
cally independent of the quark mass, while for growing n the dependence on the quark 
mass increases. As the lower moments of the distribution functions are dominated by 
the small-x region, this means that at small x quark mass effects are negligible. At 
intermediate and large x, on the other hand, the distribution functions depend strongly 
on the magnitude of the quark mass. In the limit of large quark masses the moments 
approach the predictions of the non-relativistic quark model. In particular we find 
(^x^-i^u) ^ 2(a;'^-i)('^) for all n. 

In Fig. 4 we have compared our results with the phenomenological valence quark 
distribution functions [10]. For the lowest moment the lattice result turns out to be 
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significantly larger than the phenomenological value, while for the largest moment we 
find the opposite situation. This holds for both, u and d quark distributions. Thus the 
lattice calculation predicts a valence quark distribution that is larger at intermediate 
to small X and smaller at large x than the phenomenological distribution function. At 
present we have no explanation for this discrepancy. 

Let us now discuss Aq. Sea quark effects may be neglected for heavy quarks, and 
they drop out in the difference Am — Ad. In the chiral limit we obtain 



This is to be compared with the experimental value of the axial vector coupling constant 
gA = 1.26. A recent lattice calculation of the sea quark contribution [12] finds Au = 
Ad = —0.14(5), As = —0.13(4), where we have used the perturbative renormalization 
factor. If we add these numbers to our results, we obtain for the total quark spin 
contribution to the nucleon spin 



This is in agreement with the result of a full QCD calculation [13] which includes 
dynamical quark loops. 

4 Conclusion 

We have presented some results of a high statistics calculation of the lower moments 
of the polarized and unpolarized deep-inelastic structure functions of the nucleon. The 
calculation has been performed in the quenched approximation, and it was done for 
three different quark masses. This allowed us to extrapolate our results to the chiral 
limit. 

The lattice data are rather accurate now, so that it is equally important to deter- 
mine the renormalization constants precisely. We have computed the renormalization 
constants in perturbation theory to one loop order as well as non-perturbatively. So 
far we find consistent results. 

The valence quark distributions that we have obtained on the lattice differ from 
the phenomenological distribution functions. This was also observed before [14]. One 
explanation could be that at smaller values of higher twist contributions are non- 
negligible, which have not been included in the phenomenological analysis. We plan 
to investigate this possibility in the future. Our results for Aq are consistent with 
experiment. 

It is interesting to see how the results vary with the quark mass. At large quark 
masses our results agree largely with what one would expect on the basis of the quark 
model. For small quark masses there are, however, significant changes. 



Au - Ad = gA = 1.07(9). 
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